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^ ABSTRACT 

' We show how to improve constraints on Q,m, os, and the dark-energy equation-of-state parameter, w, obtained by Mantz et al. 
; (2008) from measurements of the X-ray luminosity function of galaxy clusters, namely MACS, the local BCS and the REFLEX 
galaxy cluster samples with luminosities L > 3 x 10^^ erg/s in the 0.1-2.4 keV band. To this aim, we use Tinker et al. (2008) 
^ , mass function instead of Jenkins et al. (2001) and the M-L relationship obtained from Del Popolo (2002) and Del Popolo et 
$_( ■ al. (2005). Using the same methods and priors of Mantz et al. (2008), we find, for a ACDM universe, Qm ~ 0.28to o4 and 
; o-g = 0.78+0 05 '^l^ile the result of Mantz et al. (2008) gives less tight constraints Qm = 0.28lo:o7 and as = 0.78+o.}3- In the 
C3 , case of a wCDM model, we find Qm = 0.27to o6, as = 0.8lto.o6 a^^d w = -1.3j:o;4, while in Mantz et al. (2008) they are again 
less tight flm = 0.24^0 ^7, as = 0.85to 20 and w = -1.4lo.7. Combining the XLF analysis with the /ja^-l-CMB-l-SNIa data 
set results in the constraint Qm = 0.269 ± 0.012, erg — 0.81 ± 0.021 and w — —1.02 ± 0.04, to be compared with Mantz et al. 
'k^ . (2008), Qm = 0.269 ± 0.016, erg = 0.82 ± 0.03 and w = -1.02 ± 0.06. The tightness of the last constraints obtained by Mantz 
' et al. (2008), are fundamentally due to the tightness of the /gas +CMB-|-SNIa constraints and not to their XLF analysis. Our 
findings, consistent with w — —1, lend additional support to the cosmological-constant model. 
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1. Introduction ues of (Henry & Arnaud 1991; Reiprich & Boringer 

2002; Schuecker et al. 2003) -a result since then con- 
Cluster of galaxies are the largest gravitationally-collapsed ^^.j^ed by cosmic microwave background (CMB) studies, 
structures in the Universe. Even at the present epoch they cosmic shear, and other experiments (Spergel et al. 2007; 
are relatively rare, with only a few percent of galaxies be- Komatsu et al. 2008; Dunkley et al. 2008; Benjamin et 
ing in clusters. In the hierarchical collapse scenario for 2007; Fu et al. 2008). For precision's sake, cluster sur- 
structure formation in the universe, the number density ^^yg ^^le local universe are particularly useful for con- 
of collapsed objects as a function of mass and cosmic straining a combination of the matter density parame- 
time is a sensitive probe of cosmology. The galaxy clusters ^er 17„ and the normalization of the power spectrum of 
that occupy the high-mass tail of this population provide density fluctuations. Following the evolution of the clus- 
a powerful and relatively clean tool for cosmology, since gp^ce density over a large redshift baseline, one can 
their growth is predominantly determined by linear grav- ^reak the degeneracy between as and n„, (Rosati et al. 
itational processes. 2002). Recently, X-ray studies (Vikhlinin et al. 2009b) of 
Starting in the 1990's, analysis of massive clusters have the evolution of the cluster mass function at z = 0-0.8 
consistently indicated low values of flm (both from the have convincingly demonstrated that the growth of cos- 
baryonic fraction arguments (White et al. 1993) and mea- mic structure has slowed down at z < 1 due to the effects 
surements of the evolution in the cluster number den- of dark energy, and these measurements have been used 
sity (Eke et al. 1998; Borgani et al. 2001) and low val- 

^ as is the amplitude of the mass density fluctuation power 
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to improve the determination of the equation of state pa- 
rameter. Although the quoted cosmological test is very 
powerful, there are two main problems in practical appli- 
cations: first, theoretical predictions provide the number 
density of clusters of a given mass, while the mass itself 
is never the directly observed quantity. Second, a cluster 
sample is needed that spans a large-z baseline and is based 
on model-independent selection criteria]! 

Determining the evolution of the space density of clus- 
ters requires counting the number of clusters of a given 
mass per unit volume at different redshifts. Therefore, 
three essential tools are required for its application as a 
cosmological test: (a) an efficient method to find clusters 
over a wide redshift range, (b) an observable estimator of 
the cluster mass, and (c) a method to compute the se- 
lection function or equivalently the survey volume within 
which clusters are found. Observations of clusters in the 
X-ray band provide an efficient and physically motivated 
method of identification, which fulfills the three require- 
ments above. The X-ray luminosity, provides a very ef- 
ficient method for identifying clusters down to a given 
X-ray flux limit and hence within a known survey vol- 
ume for each luminosity Lj,, which uniquely specifies the 
cluster selection, is also a good probe of the depth of the 
cluster gravitational potential. For these reasons most of 
the cosmological studies based on clusters have used X- 
ray-selected samples. 

According to the three points quoted above, the 
recipe for constraining cosmological parameters by 
means of clusters is composed of three ingredients: 1) 
The predicted mass function of clusters, n{M,z), as a 
function of cosmological parameters {as, ^m, w, etc.). 
2) Sky surveys with well understood selection functions 
to find clusters, as well as a relation linking cluster mass 
with an observable. A successful solution to the former 
requirement has been to identify clusters by the X-ray 
emission produced by hot intracluster gas, notably using 
data from ROSAtH. 3) A tight, well-determined scaling 
relation between survey observable (e.g. L^) and mass, 
with minimal intrinsic scatter. 

Early attempts to use evolution of the cluster mass 
function as a cosmological probe were limited by small 
sample sizes and either poor proxies for the cluster mass 

^ This is so that the search volume and the number density 
associated with each cluster are uniquely identified. 

3 The ROSAT Brightest Cluster Sample (BCS; Ebeling et 
al. 1998, 2000) and ROSAT-ESO Flux Limited X-ray sample 
(REFLEX; Bohringer et al. 2004) together cover approximately 
two-thirds of the sky out to redshift z ~ 0.3 and contain more 
than 750 clusters. The Massive Cluster Survey (MACS; Ebeling 
et al. 2001, 2007) extends these data to z ~ 0.7. The ROSAT 
160 sq. degree survey, described for the first time by Vikhlinin 
et al. (1998) is a serendipitous cluster catalogue containing 
201 groups/clusters, while the ROSAT 400 sq. degree survey 
is based on 1610 high Galactic latitude ROSAT PSPC point- 
ings (Burenin et al. 2007) and includes 266 optically confirmed 
galaxy clusters, groups and individual elliptical galaxies. 



(e.g., the total X-ray flux) or inaccurate measurements 
(e.g. temperatures with large uncertainties) 

Until some years ago the obtained results for fi,n were 
several times in disagreement. Study by different authors 
(Bahcall, Fan & Gen (1997), Bahcah & Fan (1998), Sadat, 
Blanchard & Oukbir (1998), Blanchard, Bartlett & Sadat 
(1998), Blanchard & Bartlett (1998), Eke et al. (1998), 
Viana & Liddle (1999), Reichart et al. (1999), Donahue & 
Voit 1999, Borgani et al. 2001) found values for il„i span- 
ning the entire range of acceptable values: 0.2 < H.^ < 1 
(see Reichart et al. 1999). It is interesting to note that 
Viana & Liddle (1999) using the same data set as Eke et 
al. (1998) showed that uncertainties both in fitting local 
data and in the theoretical modeling could significantly 
change the final results: they found fim — 0.75 as a pre- 
ferred value with a critical density model acceptable at 
< 90% c.l. while Eke et al. (1998) found = 0.45 ± 0.2. 

The reasons leading to the quoted discrepancies have 
been studied in several papers (Eke et al. 1998; Reichart 
et al. 1999; Donahue & Voit 1999; Borgani et al. 2001) and 
can be summarized as due to: 1) The inadequate approx- 
imation given by the mass function used (e.g., Bryan & 
Norman 1998). 2) Inadequacy in the structure formation 
as described by the spherical model leading to changes in 
the threshold parameter 6c (e.g., Governato et al. 1999). 3) 
Inadequacy in the M-T relation obtained from the virial 
theorem (see Voit & Donahue 1998; Del Popolo 2002). 
4) Effects of cooling flows. 5) Determination of the X- 
ray cluster catalog's selection function. 5) Missing high 
redshift clusters in the data used (e.g., the EMSS). 6) 
Evolution of the L-T relation (Voit & Donahue 1998). 
7) The use of different best fitting procedures to get the 
constraints (Eke et al. 1998). 8) Other effects described 
in more recent papers (e.g., Mantz et al. 2008 (hereafter 
M08); Vikhlinin et al 2009b). 

The situation with the cluster mass function data has 
been dramatically improved in the past years. A large sam- 
ple of sufficiently massive clusters extending to z ~ 0.9 
has been derived from ROSAT PSPC pointed data cover- 
ing 400 deg^ (Burenin et al. 2007). Distant clusters from 
the 400d sample were then observed with Ghandra, pro- 
viding high-quality X-ray data and much more accurate 
total mass indicators (see Vikhlinin et al. 2009b) . Ghandra 
coverage has also become available for a complete sample 
of low-z clusters originally derived from the ROSAT All- 
Sky Survey (see Vikhlinin et al. 2009b). Results from deep 
Ghandra pointings to a number of low-z clusters have sig- 
nificantly improved our knowledge of the outer cluster re- 
gions and provided a much more reliable calibration of the 
Mtot vs. proxy relations than what was possible before. 
On the theoretical side, improved numerical simulations 
resulted in better understanding of measurement biases 
in the X-ray data analysis (Nagai et al. 2007; Rasia et al. 
2006; Jeltema et al. 2007). 

In the present paper, we want to show how tighter 
constraints can be obtained in M08 model improving the 
mass function adopted by them, and the scaling laws used 
(e.g., the M-T and M-L relationships). In this paper, we 
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use the observed X-ray luminosity function to investigate 
two cosmological scenarios, assiiming a spatially flat met- 
ric in both cases: the first includes dark energy in the 
form of a cosmological constant (ACDM); the second has 
dark energy with a constant cquation-of-statc parameter, 
w (wCDM). The theoretical background for this work is 
reviewed in Section 2. Section 3 presents the results and 
Section 4 the conclusions. 

2. Theory 

In the introduction, we discussed the ingredients needed 
in the recipe used to constrain cosmological parameters 
from X-ray observations. In this section, I derive an ex- 
pression for the X-Ray luminosity function (XLF) (using 
now the mass function obtained in Del Popolo (2006a, 
b)) and M-T, L-T relations obtained in Del Popolo 2002, 
and Del Popolo et al. 2005, respectively) and then I set 
some constraints to Clm, <^8 and the dark-energy equation- 
of-state parameter w, by using the data (clusters) used 
in M08, namely MACS (Massive Cluster Survey), BCS 
(Brightest Cluster Sample), and REFLEX (ROSAT ESQ 
FLUX LIMITED X-Ray SAMPLE). Following M08, the 
constraints are obtained from measurements of the X- 
Ray luminosity function of the quoted samples. The most 
straightforward mass-observable relation to complement 
these X-ray flux-limited surveys is the mass- X-ray lumi- 
nosity relation. For sufficiently massive (hot) objects at 
the relevant redshifts, the conversion from X-ray flux to 
luminosity is approximately independent of temperature, 
in which case the luminosities can be estimated directly 
from the survey flux and the selection function is identi- 
cal to the requirement of detection. A disadvantage is that 
there is a large scatter in cluster luminosities at fixed mass; 
however, sufficient data allow this scatter to be quanti- 
fied empirically. More recently, a dramatic reduction in 
luminosity-mass scatter has been demonstrated when lu- 
minosities are measured excluding cluster centers (typi- 
cally r < 0.15r5oo; Maughan 2007; Zhang et al. 2007). 
Alternative approaches use cluster temperature (Henry 
2000; Scljak 2002; Pierpaoh et al. 2003; Henry 2004), gas 
fraction (Voevodkin & Vikhlinin 2004) or Yx parameter 
(Kravtsov et al. 2006) to achieve tighter mass-observable 
relations at the expense of reducing the size of the samples 
available for analysis. The need to quantify the selection 
function in terms of both X-ray flux and a second observ- 
able additionally complicates these efforts. 

The first ingredient of the quoted recipe (i.e., mass 
function), used in M08 was the Jenkins et al. (2001) (here- 
after JOl) mass function. JOl wrote the mass function of 
galaxy clusters of mass M at redshift z as a "universal 
function" of ^-^(M,^) 



_ M n{M,z) 



which was fitted by 
/(a-i)=Ae(-l'"-"'+^l') 



(1) 



(2) 



for cosmological -constant models, with A = 0.316, B = 
0.67, and e = 3.82. 

As shown in Del Popolo (2006a) and Del Popolo 
(2006b), the theoretical mass function obtained in the 
quoted papers is in better agreement with high resolu- 
tion N-body simulations, namely Reed et al. 2003 (R03), 
Yahagi et al. (2004) (YNY), Warren et al. 2006 (W06), 
and Tinker et al. (2008) (see the following and Fig. lb) 
(T08). 

The mass function was calculated according to the 
model of Del Popolo (2006a, b) . The multiplicity function, 
in the quoted model, is given by: 



where 




{avf 



[av) 



{avf 



(3) 



g{ai) = I 1 - a, + 



ai{ai - 1) 
2! 



ai{ai - 1) • ■ • (oii - 4) 
5! 



(4) 



where i = 1 or 2, ai = 0.585, /3i = 0.46, a2 = 0.5 and 
132 = 0.35, as = 0.4 and = 0.02, a = 0.707, and Ai = 
1.2 is the normalization constant. 

The "multiplicity function" is correlated with the 
usual, more straightforwardly used, "mass function" as 
follows. Following Sheth & Tormen (2002) (hereafter ST) 
notation, if f(M, S)dM denotes the fraction of mass that is 
contained in collapsed haloes that have mass in the range 
M-M + dM, at redshift z, and S{z) is the redshift de- 
pendent overdensity, the associated "unconditional" mass 
function is: 



n{M,S)dM = ^f{M,S)dM 



(5) 



In Fig. la, we plot the multiplicity function obtained 
in this paper (simbols are described in the figure caption). 

There are some differences between the quoted sim- 
ulations and the JOl simulations. First, the multiplicity 
function of the present paper, similar to that of YNY, 
in the low-z^ region of v < 1 systematically falls below 
the JOl functions. In this region the multiplicity func- 
tion of the present paper is very close to that of YNY. 
Additionally, the numerical multiplicity functions (and 
that in Del Popolo 2006a, b) have an apparent peak at 
i/ ~ instead of the plateau that is seen in the JOl function. 
Similar differences are seen in the high-i/ region. These dif- 
ferences between mimerical multiplicity functions (R03; 
YNY; W05; Del Popolo 2006a,b) and JOl, are however 
within 1-cr error-bars, and so they are overall in agree- 
ment. The multiplicity function obtained in the present 
paper has a peak at i/ ~ 1 as in YNY numerical multiplic- 
ity function, instead of a plateau as in the JOl function. 
Differences are observed also in the redshift evolution of 
the JOl mass function (Del Popolo 2006b). Summarizing 
the fitting formulas presented by JOl are accurate to 
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2± 10 - 20% (Tinker ct al. 2008 (T08)). In our model, the 
mass function that we used is given by Eq. Eq.(|4|), 
and Eq© which is in perfect agreement with the T08 
mass function, as shown in Fig. (lb). So, the accuracy of 
the mass function is, as in T08, of the order of ~ 5% for 
ACDM models for the mass and redshift range of interest 
in this study. As a consequence, in this way the theoret- 
ical uncertainties in the mass function do not contribute 
significantly to the systematic error budget. 

In Fig. lb, I plot the mass function for all of our out- 
puts in the /(cr) — \n{a^^) plane. Large values of Incr^^ 
correspond to rare haloes of high redshift and/or high 
mass, while small values of lncr~^ describe haloes of low 
mass and redshift combinations. Fig. lb shows the func- 
tion /(cr) measured for all simulations in Table 1 of T08, 
the solid line the fit to the data (namely T08 eq. 3) and 
the dashed line the model of the present paper. 

As previously reported, one of the main problems of us- 
ing the mass function to constrain cosmological parameter 
is that theoretical predictions provide the number density 
of clusters of a given mass, while the mass itself is never 
the directly observed quantity. One then needs relations 
connecting mass with other quantities more easily obtain- 
able which can be used as a surrogate for cluster mass. 
Over the past decade, observations of clusters of galax- 
ies (e.g. ROSAT, ASCA) have shown the existence of a 
correlation between the total gravitating mass of clusters, 

O 

their X-ray luminosity (Lx) and the temperature 
(Tx) of the intracluster medium (ICM) (David et al. 1993; 
Markevitch 1998; Horner, Mushotzky & Scharf 1999). By 
means of the quoted scaling relations one can obtain differ- 
ent methods for tracing the evolution of the cluster num- 
ber density: (1) The X-ray temperature function (XTF), 
which has been presented for local (e.g., Henry & Arnaud 
1991) and distant clusters (Eke et al. 1998; Henry 2000). 
2) The evolution of the X-ray luminosity function (XLF). 
In this case, we need a relation between the observed Lx 
and the cluster virial mass. 

In the following, following M08, we shall use the XLF 
to constrain cosmological parameters. Then the next cru- 
cial step, after having a mass function, is to convert it 
in a Luminosity function (XLF). This can be done by 
first converting mass into intra-cluster gas temperature, 
by means of the M — relation, and then converting 
the temperature into X-ray luminosity, by means of the 
Lx — Tx relation. MOB used a self-similar relationship be- 
tween mass and X-ray luminosity for massive clusters 
(e.g., Bryan & Norman 1998) modified by an additional 
redshift-dependent factor (see Morandi et al. 2007). At 
this point, we must stress an important point. Numerical 
simulations confirm that the DM component in clusters 
of galaxies, which represents the dominant fraction of the 
mass, has a remarkably self-similar behavior; however the 
baryonic component does not show the same level of self- 

* Since Mtot compares with the ICM temperature measure- 
ments that can be obtained through X-ray spectroscopy, this 
explains the importance of a mass-temperature (MT) relation. 



similarity. This picture is confirmed by X-ray observa- 
tions, see for instance the deviation of the L-T relation 
in clusters, which is steeper than the theoretical value 
predicted by the previous scenario. More precisely, until 
some years ago, the cluster structure was considered to 
be scale-free, which means that the global properties of 
clusters, such as halo mass, luminosity-temperature, and 
X-ray luminosity would scale self-similarly (Kaiser 1986; 
Evrard & Henry 1991). In particular, the gas temperature 
would scale with cluster mass as T esc M^/^ and the bolo- 
metric X-ray luminosity would scale with temperature as 
i oc , in the bremsstrahlung-dominated regime above 2 
keV. Studies following that of Kaiser (1986) showed that 
the observed luminosity-temperature relation is closer to 
L (xT^ (e.g., Edge & Stewart 1991), indicating that non- 
gravitational processes should influence the density struc- 
ture of a clusters core, where most of the luminosity is 
generated (Kaiser 1991; Evrard & Henry 1991; Navarro 
et al. 1995; Bryan & Norman 1998). One way to obtain a 
scaling law closer to the observational one is to have non- 
gravitational energy injected into the ICM before or dur- 
ing cluster formation, the so-called pre-heating (Ponman 
et al. 1999; Bower et al. 1997; Cavaliere et al. 1997, 1999; 
Tozzi & Norman 2001; Borgani et al. 2001; Voit & Brian 
2001), feedback processes that alter the gas characteristics 
during the evolution of the cluster (Voit & Bryan 2001), 
cooling flows (Allen & Fabian 1998). A similar situation 
is valid for the M-T relationship, namely that the self- 
similarity in the M-T relation seems to break at a few 
keV (Nevalanien et al. 2000; Xu, Finoguenov et al. 2001; 
Muanwong et al. 2001; Bialek, Evrard & Mohr 2001). 
Consequently, if, as in M08, one starts with self-similar 
scaling laws in order to have consistent scaling relations 
one has to compare the self-similar scaling relations to 
observations (Morandi et al. 2007). 

Different from the M08 approach, in the following, we 
use models for the L-T, T-M, relationships taking into 
account the non-self similarity: namely, the M — rela- 
tion obtained analytically using the model of Del Popolo 
(2002), while the — relation is that obtained in Del 
Popolo, Hiotelis & Peharrubia (2005) based on an im- 
provement of the Punctuated Equilibrium Model (PEM) 
of Cavaliere et al. (1997, 1998, 1999). The drawbacks of 
using self-similar relationships fitted to the data (clusters) 
and the reasons to use a different approach, were already 
discussed in Del Popolo (2003) (their sect. 3), and in the 
remainder of this section. 

Similarly, to the present study, in Del Popolo (2003) 
we used the models for the L-T, T-M, relationships 
instead of the scaling relations obtained from simulations 
of Chandra data (see, e.g., Pierpaoli et al. 2001, 2003 
for references) H Eq. (5) in M08 similar to that Eq. 
(13) of Pierpaoli et al. (2001) or Eq. (4) of Pierpaoh 
et al. (2003) comes from rather simplistic arguments 
(dimensional analysis and an assumption that clusters are 



^ Notice that Eq. 22 in Del Popolo 2003, and Eq. 5 in M08, 
are very similar. 
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self-similar) not taking into account important physical 
effects that gives rise to a non-self-similar behavior of 
the quoted relation, as previously discussed. The fitting 
procedure used by M08, trying to take account of the 
previous physics and the non-sclf-similar behavior of the 
relationship, is complicated by several effects. In fact, 
in the fit one uses data that may contain small groups 
which can be influential in the estimation of the slope 
of the model, and one has then to choose accurately 
the data to be used in the fit. This choice mitigates 
the possibility of obtaining biased results if slope of the 
mass-luminosity relation is different for massive clusters 
compared with smaller groups. In M08 they fitted only 
the data (clusters) with L > 3 x 10'*'* crg/s in the 
0.1-2.4 keV band. Moreover, the process of fitting the 
model in Eq. 7 of M08 is complicated by the presence 
of Malmquist bias. Close to the flux limit for selection, 
any X-ray selected sample will preferentially include the 
most luminous sources for a given mass. This results 
in a steepening of the derived mass-luminosity relation 
and a bias in the inferred intrinsic scatter in luminosity 
for a given mass. The use of the extended sample of 
Rciprich & Boringcr (2002) (RB02), rather than only 
their flux-limited HIFLUGCS sample, partially mitigates 
this effect by softening the flux limit. A further problem 
is that as a consequence of Malmquist bias there is a 
strong apparent, but not necessarily physical, correlation 
between luminosity and redshift due to the fact that the 
flux limit corresponds to higher luminosities at higher 
redshifts. 



For what concerns the data (clusters) used in the anal- 
ysis, they arc the same of those used by MOB: the follow- 
ing three flux-limited surveys are included in our analy- 
sis: the BCS (EbeUng et al. 1998) and REFLEX sample 
(Bohringer et al. 2004) at low redshifts {z < 0.3), and the 
MACS (Ebeling et al. 2001) at 0.3 < z < 0.5 (see M08). In 
the analysis, the sample was chosen to cover the redshift 
range z < 0.5, since at higher redshifts the number of un- 
relaxed clusters decrease, and the L-T and T-M relations 
are appropriate for relaxed clusters. The purpose of this 
paper is to present an analysis based only on the X-ray lu- 
minosity function (XLF) data described above, along with 
the priors described in Sect. 4 of M08. 

Following M08, we parametrize the full model flt- 
ted to the X-ray luminosity function data as h, ^5/1^, 
rich^, as, rig, w where flj, and ilg are the baryon and 
cold dark matter densities {Sim = fib + Sic)- In addi- 
tion to the assumption of spatial flatness, we adopt the 
Gaussian priors h = 0.72 ± 0.08 (Freedman et al. 2001) 
and Qbh^ = 0.0214± 0.002 (Kirkman et al. 2003) from the 
Hubble Key Project and Big Bang nucleosynthesis stud- 
ies, respectively. Since the results are insensitive to the 
spectral index within a reasonable range (see M08), we flx 
Hs = 0.95 in accordance with (Spcrgcl et al. 2007) for the 
standard analysis. The dark-energy equation of state was 
bounded by a uniform prior, — 5 < w < 0. 



The luminosity function likelihood is the same as in 
MOB (Sect. 4.2). 

The likelihood that N clusters with inferred luminosi- 
ties in a range dL exist in a volume dV can in general 
be written as a Poisson probability plus a correction due 
to the clustering of halos with one another. If the plane 
of redshift and inferred luminosity is divided into non- 
overlapping cells, then the likelihood of our data is simply 



P{N,,N2,...) = ll 



(6) 



where Nj and Nj are the number of clusters detected and 
predicted in the jth cell, respectively. 

If the cells are taken to be rectangular, with the j'th 
cell given by z-^^ < z < zf^ and L^^^ < L < L'f \ then 



(2) f(2) 

dz Jtw 



dn{z, L) 
dL 



(7) 



where V{z) is the comoving volume within redshift z. In 
the absence of intrinsic scatter in the mass-luminosity re- 
lation and measurement errors in the observed luminosi- 
ties, the derivative of the comoving number density would 
be simply 



dh{z, L) 
dL 



fsky{z, L) 



dM{L)dn{z,M) 



dL 



dM 



(8) 



Here fsky is the sky coverage fraction of the surveys as 
a function of redshift and inferred luminosity, dn/dM is 
no longer the Jenkins mass function but the one disciissed 
in the present paper and M{L) is the mass-luminosity 
relation discussed in the present paper. 

Similar to MOB, the presence of scatter requires us to 
take into account that a cluster detected with inferred lu- 
minosity L could potentially have any true luminosity L 
and mass M, with some associated probability. To cal- 
culate the predicted number density correctly, we must 
therefore convolve with these probability distributions: 



dn{z, L) 
dL 



— fsky {z 



f 

Jo 



,L) / 
Jo 



dL P{L\L) 



(9) 



dM P{L\M) 



dn{z, M) 
dM ■ 



P{L\M) is a log-normal distribution whose width is like in 
MOB the intrinsic scatter in the mass-luminosity relation, 
r}{z), and P{L\L) is a normal distribution whose width as 
a function of flux is modeled as a power law, as described 
in Sect. 3.2 of MOB. 



3. Results 

In Fig. 2, we compare, for a ACDM the joint Slm,-a% con- 
straints obtained from the BCS, REFLEX and MACS 
data sets combination The marginalized constraints from 
the combination of the three cluster samples are Sim = 
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0.281^;^^ and erg = 0-78 om while the result of M08 gives 
less tight constraints flm — 0.28]to o7 a-nd erg ~ 0.78lg ijg. 

Our previous constraints are in good agreement with 
recent, independent results from the CMB (Spergel et al. 
2007) and cosmic shear, as measured in the 100 Square 
Degree Survey (Benjamin et al. 2007) and CFHTLS Wide 
field (Fu et al. 2008). Our results are also in good overall 
agreement with previous findings based on the observed X- 
ray luminosity and temperature functions of clusters (Eke 
et al 1998, Donahue & Voit 1999, Henry 2000, Borgani 
et al. 2001, Seljak 2002, Allen et al. 2003, Pierpaoh et 
al. 2003, Schuecker et al. 2003, Henry 2004). Our resuh 
on is in excellent agreement with current constraints 
based on cluster fgas data (Allen et al. 2008 and references 
therein) and the power spectrum of galaxies in the 2dF 
galaxy redshift survey (Cole et al. 2005) and Sloan Digital 
Sky Survey (SDSS) (Eisenstein et al. 2005, Tegmark et al. 
2006, Percival et al. 2007), as well as the combination of 
CMB data with a variety of external constraints (Spergel 
2007). 

In Fig. 3a, we set constraints for the wCDM model, and 
we plot the joint constraints on Vim and erg from the lumi- 
nosity function data using our standard priors, while Fig. 
3b displays constraints on i7,„ and w obtained indepen- 
dently from the XLF data. The marginalized results from 
the X-ray luminosity function data are flm — 0.27^Qgg, 
CTg = 0.8llJ^{!^ and w = -1.3t^j, while in M08 they 
are again less tight r2„ = 0.24to:o7, o-g = 0.85lo;2o and 
w = — 1.4^q'7. Our new XLF results are consistent with 
the cosmological-constant model {w = —1). 

An improvement on the previous results can 
be obtained by adding to the XLF analysis the 
/gas+CMB-f SNIa data set. The /gas+CMB-hSNIa combi- 
nation already provides tight constraints on ftm, h, Vlbh? 
and Ug (hence no priors on these parameters are used in 
either combined analysis), but the degeneracy between w 
and (Tg (right panel of Fig. 4) limits the precision of the 
dark energy results. The addition of the XLF data breaks 
the degeneracy in the ilm-crg plane (left panel), resulting 
in tighter constraints on fim, fs and w. The degeneracy 
breaking power of other combinations of data with the 
CMB is discussed by Spergel et al. (2007). The result- 
ing constraint are Vt^ = 0.269 ± 0.012, erg ^ 0.81±0.021 
and w — — 1.02±0.04, to be compared with M08 Q.m = 
0.269±0.016, CTg = 0.82±0.03 and w = -1.02±0.06. The 
previous constraints are in agreement within the errors 
with Vikhlinin et al. (2009b) constraints, namely ftm — 
0.255 ± 0.043, CTg = 0.786±0.011 and w = -0.991±0.045. 

It is important to note that the tight constraints ob- 
tained by M08 when combining XLF analysis with the 
/gas+CMB+SNIa data set are primarily due to the tight- 
ness of the constraints obtained from /gag-l-CMB-l-SNIa 
data itself and not to the precision of the XLF analysis 
of M08, as shown by comparing our results for flm, CTg, w 
obtained by the XLF analysis with those of M08. In our 
model the improvement in the mass function model and 



the L-M relationship gives rise to tight constraints even 
when using only the XLF function. 

In order to understand why the results of our analysis 
are different from those of M08, we have to stress a key 
point. The M08 paper, as well as several others papers 
in the literature, used two different data sets: REFLEX, 
BCS, and MACS to constrain cosmological parameters 
and another external data set to constrain the luminosity- 
mass relation (RB02 data set), which in M08 is a power- 
law with three free parameters, without explicitly account- 
ing for selection bias. Consequently, it was necessary to re- 
strict that external data set (RB02) to low redshifts and 
high fluxes in order to minimize the effects of selection 
bias, making it impossible to test for departures from self- 
similar evolution in the scaling relation. In order to have a 
"self-consistent" analysis, it is necessary that a single like- 
lihood function be applied to the full data set which en- 
compasses the entire theoretical model (cosmology -I- scal- 
ing relations) so as to ensure that the covariance among 
all the model parameters is fully captured and that the 
effects of the mass function and selection biases are prop- 
erly accounted for throughout. This kind of analysis was 
performed for the first time by Vikhlinin et al. (2009a, b), 
who used the same cluster sample to constrain the scaling 
relations, thus obtaining tighter constraints. 

In the analysis of the present paper, the L-M relation 
is a physically motivated relation (not a power-law with 
free parameters) which does not require fits to data, as 
in M08. Since we do not need the double analysis of M08 
and previous papers, the first to get the L-M fitting pa- 
rameters from RB02 data, and the second to obtain the 
cosmological constraints using BCS, REFLEX and MACS, 
we bypass the quoted drawback in the M08 analysis. 

It is interesting to note that a month after the present 
paper was submitted, two papers, Mantz et al. (2009a, b), 
appeared in arXiv showing that the key point that I pre- 
viously stressed, namely generalizing M08 to allow the 
quoted simultaneous and self-consistent fit and using T08 
mass function (instead of that in Jenkins et al. 2001) re- 
sult in cosmological constraints that are a factor 2-3 better 
than those in M08, based on the same flux-limited sample 
of clusters. In the present paper, we have also checked that 
using the same L-M relation used in M08, we reobtain the 
same set of constraints derived by MO^ . 

Another point to stress concerns the use of our non- 
self similar L-M relation for clusters of luminosity L > 
3 X 10^'* erg/s in the 0.1-2.4 keV band. Since the clusters 
included in the M08 sample are high X-ray luminosity 
(above 3 x 10^^ erg/s), one could think that the changes 
in the L-M relation of the present paper, with respect to 
the classical self-similar model, will not produce signifi- 
cant changes in constraints on the cosmological parame- 
ters. Even if major differences between the L-M model 
of the present paper and the self-similar model are ob- 
served at gas temperatures below 3 keV, we stress that 

" Mantz et al. (2009a,b) obtain r2„ = 0.27 ± 0.02, erg = 
0.79±0.03 and w = -0.96±0.06. 
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the present L-M relation depends on the M-T and L-T 
relationships, and especially the second one (based on the 
Modified Punctuated Equilibrium Model (MPEM)) never 
behaves in a self-similar way as shown in Del Popolo et 
al. (2005) (even at gas temperatures higher than 10 keV). 
Moreover, as previously reported, the improvement in the 
constraints is strictly connected to the fact that we bypass 
the quoted drawback in the M08 analysis by means of our 
L-M relation not depending on parameters that must be 
fixed using external data. 

In order to obtain tighter and tighter constraints one 
needs to try to reduce to the minimum the systematic un- 
certainties in the analysis. Much progress is expected over 
the coming years in refining the ranges of these allowances, 
both observationally and through improved simulations. A 
reduction in the size of the required systematic allowances 
will tighten the cosmological constraints. Improved nu- 
merical simulations of large samples of massive clusters, 
including a more complete treatment of star formation and 
feedback physics that reproduces both the observed opti- 
cal galaxy luminosity function and cluster X-ray proper- 
ties, will be of major importance. Further deep X-ray and 
optical observations of nearby clusters will provide better 
constraints on the viscosity of the cluster gas. Improved 
optical/near infrared observations of clusters should pin 
down the stellar mass fraction in galaxy clusters and its 
evolution. Ground and space-based gravitational lensing 
studies will provide important, independent constraints on 
the mass distributions in clusters; a large program using 
the Subaru telescope and Hubble Space Telescope is un- 
derway, as is similar work by other groups (e.g. Hoekstra 
2007). 

In the near future, continuing programs of Chandra 
and XMM-Newton observations of known. X-ray lumi- 
nous clusters should allow important progress to be made, 
both by expanding the fgas sample (e.g. Chandra snap- 
shot observations of the entire MACS sample; Ebeling et 
al. 2001, 2007) and through deeper observations of the 
current target list. A new, large area X-ray survey such 
as that proposed by the Spectrum- RG/eROSITA project 
could make a substantial contribution, finding hundreds 
of suitable systems at high redshifts. 

Looking a decade ahead, the International X-ray 
Observatory (IXO), result of the merging of NASA's 
Constellation-X and ESA/JAXA's XEUS mission con- 
cepts, will offer the possibility to carry out precise studies 
of dark energy using the fgas technique^ 



The large collecting area and combined spatial/spectral re- 
solving power of IXO should permit precise /gas measures. An 
investment of ~ 10 Ms of IXO time to measure fgas to 5% 
(corresponding to 3.3% accuracy in distance) in each of the 
500 hottest, most X-ray luminous, dynamically relaxed clus- 
ters detected in future cluster surveys, spanning the redshift 
range < z < 2 (typical redshift z ~ 0.6), will be sufficient 
to constrain cosmological parameters with a DETF figure of 
merit of 2040. 



4. Conclusions 

In the present paper, we showed how to improve the con- 
straints on rim, cg, and the dark-energy equation-of-state 
parameter, w, from measurements of the X-ray luminos- 
ity function of galaxy clusters, as performed by M08. 
Improving the mass function by means of Del Popolo 
(2006a, b) model, which was shown to be in good agree- 
ment with T08 and using the L-M relationship obtained in 
Del Popolo (2002) and Del Popolo et al. (2005), we showed 
that the XLF alone can give tight constraints on the cos- 
mological parameters. Using the same methods and priors 
of M08, we find, for a ACDM universe, = 0.28toM and 
(Ts = 0.78^Q Q5 and similarly in the case of a wCDM model, 
we find = 0.271° ag = OMt^o fe and w = -l/St^ol 
both tighter than M08 results. Combining the XLF anal- 
ysis with the /gas+CMB+SNIa data set results in the 
constraint = 0.269 ± 0.012, erg = 0.81 ± 0.021 and 
w = — 1.02 ± 0.04, in agreement with the most recent de- 
termination of the quoted parameters (Allen et al. 2008; 
Vikhhnin et al. 2009b; Percival et al. 2009). Our findings, 
consistent with w — —1 lends additional support to the 
cosmological-constant model. 
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Fig. 1. In Panel (a), the solid line represents the multi- 
plicity function obtained in this paper, the short-dashed 
line the fitting formula proposed by Yahagi et al. (2004) 
(their Eq. 7), the dotted line the Sheth & Tormen 
(2002) (ST) multiplicity function, the long-dashed line the 
Jenkins et al. (2001) multiplicity function. The errorbars 
with open circles represent the run 140 of YNY, those 
with filled squares the case 70b, those with open squares 
the case 70a, those with filled circles the case 35b, those 
with crosses the case 35a. Panel (b) . Mass function plotted 
in redshift-independent form. The measured /(cr) from all 
simulations in Table 1 of Tinker et al. (2008). The solid 
line is the best fit function of equation (3) (Tinker et al. 
2008). The dashed line the model in the present paper. 




Fig. 2. Joint 68.3 and 95.4 per cent confidence constraints 
on n„, and erg for a ACDM model from MACS, BCS, and 
REFLEX combination using standard priors (as described 
in the text). 
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Fig. 3. Panel (a). Same as Fig. 2 but for a wCDM model. 
Panel (b). Joint 68.3 and 95.4 per cent confidence con- 
straints on flm and w for a constant-w model using the 
X-ray luminosity function data and standard priors (as 
described in the text). 
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Fig. 4. Joint 68.3 and 95.4 per cent confidence constraints 
on rtm and as (left panel) and cts and w (right panel) ob- 
tained from a combined /gas+CMB+SNIa analysis (blue) 
and the improved constraints obtained by combining these 
data with the XLF (gold). No priors on h, n^h^ or Us 
are imposed in cither analysis. In the left panel, the re- 
sults from the XLF alone using standard priors are shown 
(purple) in order to illustrate the degeneracy breaking. 
Note that in the left panel we plotted just the inner confi- 
dence contours in the /ga^+CMB-l-SNIa analysis, in order 
to have a more readable plot. 



